DeTox: a pipeline for the detection of toxins in venomous organisms.

Venomous organisms have independently evolved the ability to produce toxins 101 times during their evolutionary history, resulting in over 200 000 venomous species. Collectively, these species produce millions of toxins, making them a valuable resource for bioprospecting and understanding the evolutionary mechanisms underlying genetic diversification. RNA-seq is the preferred method for characterizing toxin repertoires, but the analysis of the resulting data remains challenging. While early approaches relied on similarity-based mapping to known toxin databases, recent studies have highlighted the importance of structural features for toxin detection. The few existing pipelines lack an integration between these complementary approaches, and tend to be difficult to run for non-experienced users. To address these issues, we developed DeTox, a comprehensive and user-friendly tool for toxin research. It combines fast execution, parallelization and customization of parameters. DeTox was tested on published transcriptomes from gastropod mollusks, cnidarians and snakes, retrieving most putative toxins from the original articles and identifying additional peptides as potential toxins to be confirmed through manual annotation and eventually proteomic analysis. By integrating a structure-based search with similarity-based approaches, DeTox allows the comprehensive characterization of toxin repertoire in poorly-known taxa. The effect of the taxonomic bias in existing databases is minimized in DeTox, as mirrored in the detection of unique and divergent toxins that would have been overlooked by similarity-based methods. DeTox streamlines toxin annotation, providing a valuable tool for efficient identification of venom components that will enhance venom research in neglected taxa.

Relaxation of Natural Selection in the Evolution of the Giant Lungfish Genomes.

Nonadaptive hypotheses on the evolution of eukaryotic genome size predict an expansion when the process of purifying selection becomes weak. Accordingly, species with huge genomes, such as lungfish, are expected to show a genome-wide relaxation signature of selection compared with other organisms. However, few studies have empirically tested this prediction using genomic data in a comparative framework. Here, we show that 1) the newly assembled transcriptome of the Australian lungfish, Neoceratodus forsteri, is characterized by an excess of pervasive transcription, or transcriptional leakage, possibly due to suboptimal transcriptional control, and 2) a significant relaxation signature in coding genes in lungfish species compared with other vertebrates. Based on these observations, we propose that the largest known animal genomes evolved in a nearly neutral scenario where genome expansion is less efficiently constrained.

Increasing temperature-driven changes in life history traits and gene expression of an Antarctic tardigrade species.

The Antarctic region has been experiencing some of the planet's strongest climatic changes, including an expected increase of the land temperature. The potential effects of this warming trend will lead ecosystems to a risk of losing biodiversity. Antarctic mosses and lichens host different microbial groups, micro-arthropods and meiofaunal organisms (e.g., tardigrades, rotifers). The eutardigrade Acutuncus antarcticus is considered a model animal to study the effect of increasing temperature due to global warming on Antarctic terrestrial communities. In this study, life history traits and fitness of this species are analyzed by rearing specimens at two different and increasing temperatures (5°C vs. 15°C). Moreover, the first transcriptome analysis on A. antarcticus is performed, exposing adult animals to a gradual increase of temperature (5°C, 10°C, 15°C, and 20°C) to find differentially expressed genes under short- (1 day) and long-term (15 days) heat stress. Acutuncus antarcticus specimens reared at 5°C live longer (maximum life span: 686 days), reach sexual maturity later, lay more eggs (which hatch in longer time and in lower percentage) compared with animals reared at 15°C. The fitness decreases in animals belonging to the second generation at both rearing temperatures. The short-term heat exposure leads to significant changes at transcriptomic level, with 67 differentially expressed genes. Of these, 23 upregulated genes suggest alterations of mitochondrial activity and oxido-reductive processes, and two intrinsically disordered protein genes confirm their role to cope with heat stress. The long-term exposure induces alterations limited to 14 genes, and only one annotated gene is upregulated in response to both heat stresses. The decline in transcriptomic response after a long-term exposure indicates that the changes observed in the short-term are likely due to an acclimation response. Therefore, A. antarcticus could be able to cope with increasing temperature over time, including the future conditions imposed by global climate change.

Bivalves Present the Largest and Most Diversified Repertoire of Toll-Like Receptors in the Animal Kingdom, Suggesting Broad-Spectrum Pathogen Recognition in Marine Waters.

Toll-like receptors (TLRs) are the most widespread class of membrane-bound innate immune receptors, responsible of specific pathogen recognition and production of immune effectors through the activation of intracellular signaling cascades. The repertoire of TLRs was analyzed in 85 metazoans, enriched on molluscan species, an underrepresented phylum in previous studies. Following an ancient evolutionary origin, suggested by the presence of TLR genes in Anthozoa (Cnidaria), these receptors underwent multiple independent gene family expansions, the most significant of which occurred in bivalve molluscs. Marine mussels (Mytilus spp.) had the largest TLR repertoire in the animal kingdom, with evidence of several lineage-specific expanded TLR subfamilies with different degrees of orthology conservation within bivalves. Phylogenetic analyses revealed that bivalve TLR repertoires were more diversified than their counterparts in deuterostomes or ecdysozoans. The complex evolutionary history of TLRs, characterized by lineage-specific expansions and losses, along with episodic positive selection acting on the extracellular recognition domains, suggests that functional diversification might be a leading evolutionary force. We analyzed a comprehensive transcriptomic data set from Mytilus galloprovincialis and built transcriptomic correlation clusters with the TLRs expressed in gills and in hemocytes. The implication of specific TLRs in different immune pathways was evidenced, as well as their specific modulation in response to different biotic and abiotic stimuli. We propose that, in a similar fashion to the remarkable functional specialization of vertebrate TLRs, the expansion of the TLR gene family in bivalves attends to a functional specification motivated by the biological particularities of these organisms and their living environment.

Unveiling the Impact of Gene Presence/Absence Variation in Driving Inter-Individual Sequence Diversity within the CRP-I Gene Family in Mytilus spp.

Mussels (Mytilus spp.) tolerate infections much better than other species living in the same marine coastal environment thanks to a highly efficient innate immune system, which exploits a remarkable diversification of effector molecules involved in mucosal and humoral responses. Among these, antimicrobial peptides (AMPs) are subjected to massive gene presence/absence variation (PAV), endowing each individual with a potentially unique repertoire of defense molecules. The unavailability of a chromosome-scale assembly has so far prevented a comprehensive evaluation of the genomic arrangement of AMP-encoding loci, preventing an accurate ascertainment of the orthology/paralogy relationships among sequence variants. Here, we characterized the CRP-I gene cluster in the blue mussel Mytilus edulis, which includes about 50 paralogous genes and pseudogenes, mostly packed in a small genomic region within chromosome 5. We further reported the occurrence of widespread PAV within this family in the Mytilus species complex and provided evidence that CRP-I peptides likely adopt a knottin fold. We functionally characterized the synthetic peptide sCRP-I H1, assessing the presence of biological activities consistent with other knottins, revealing that mussel CRP-I peptides are unlikely to act as antimicrobial agents or protease inhibitors, even though they may be used as defense molecules against infections from eukaryotic parasites.

Independent acquisition of short insertions at the RIR1 site in the spike N-terminal domain of the SARS-CoV-2 BA.2 lineage.

Although the major SARS-CoV-2 omicron lineages share over 30 non-synonymous substitutions in the spike glycoprotein, they show several unique mutations that were acquired after their ancestral split. One of the most intriguing mutations associated with BA.1 is the presence of the inserted tripeptide Glu-Pro-Glu within the N-terminal domain, at a site that had previously independently acquired short insertions in several other SARS-CoV-2 lineages. Although the functional implications of the small nucleotide sequences found at this insertion hotspot, named RIR1, are still unclear, we have previously hypothesized that they may play a compensatory role in counterbalancing minor fitness deficits associated with other co-occurring spike non-synonymous mutations. Here, we show that similar insertion events have independently occurred at RIR1 at least 20 times in early 2022 within the BA.2 lineage, being occasionally associated with significant community transmission. One of these omicron sublineages, characterized by a Ser-Gly-Arg insertion in position 212, has been responsible for over 4000 documented COVID-19 cases worldwide between January and July 2022, for the most part concentrated in Denmark, where it reached a national prevalence close to 4% (10% in the Nordjylland region) in mid-May. Although the concurrent spread of the BA.2.12.1, BA.4 and BA.5 lineages led to the rapid decline of this BA.2 sublineage, the independent acquisition of several other RIR1 insertions on a BA.2 genomic background suggests that these events may provide a slight fitness advantage. Therefore, they should be carefully monitored in the upcoming months in other emerging omicron-related lineages, including BA.5.

Detecting Structural Variants and Associated Gene Presence-Absence Variation Phenomena in the Genomes of Marine Organisms.

As complete genomes become easier to attain, even from previously difficult-to-sequence species, and as genomic resequencing becomes more routine, it is becoming obvious that genomic structural variation is more widespread than originally thought and plays an important role in maintaining genetic variation in populations. Structural variants (SVs) and associated gene presence-absence variation (PAV) can be important players in local adaptation, allowing the maintenance of genetic variation and taking part in other evolutionarily relevant phenomena. While recent studies have highlighted the importance of structural variation in Mollusca, the prevalence of this phenomenon in the broader context of marine organisms remains to be fully investigated.Here, we describe a straightforward and broadly applicable method for the identification of SVs in fully assembled diploid genomes, leveraging the same reads used for assembly. We also explain a gene PAV analysis protocol, which could be broadly applied to any species with a fully sequenced reference genome available. Although the strength of these approaches have been tested and proven in marine invertebrates, which tend to have high levels of heterozygosity, possibly due to their lifestyle traits, they are also applicable to other species across the tree of life, providing a ready means to begin investigations into this potentially widespread phenomena.

The Mantle Transcriptome of Chamelea gallina (Mollusca: Bivalvia) and Shell Biomineralization.

The striped venus Chamelea gallina is a bivalve mollusc that represents one of the most important fishery resources of the Adriatic Sea. In this work, we investigated for the first time the ability of this species to modulate the expression of genes encoding proteins involved in biomineralization process in response to biotic and abiotic factors. We provided the first comprehensive transcriptome from the mantle tissue of clams collected in two sampling sites located along the Italian Adriatic coast and characterized by different environmental features. Moreover, the assessment of environmental parameters, scanning electron microscopy (SEM), and X-ray diffraction (XRD) measurements on valves were conducted to better contextualize RNA sequencing (RNA-Seq) data. Functional annotation of differentially expressed genes (DEGs) and SEM observations highlighted a different shell mineralization behaviour in C. gallina clams collected from two selected sites characterized by diverse environmental parameters.

Transcriptional Contribution of Transposable Elements in Relation to Salinity Conditions in Teleosts and Silencing Mechanisms Involved.

Fish are an interesting taxon comprising species adapted to a wide range of environments. In this work, we analyzed the transcriptional contribution of transposable elements (TEs) in the gill transcriptomes of three fish species exposed to different salinity conditions. We considered the giant marbled eel Anguilla marmorata and the chum salmon Oncorhynchus keta, both diadromous, and the marine medaka Oryzias melastigma, an euryhaline organism sensu stricto. Our analyses revealed an interesting activity of TEs in the case of juvenile eels, commonly adapted to salty water, when exposed to brackish and freshwater conditions. Moreover, the expression assessment of genes involved in TE silencing mechanisms (six in heterochromatin formation, fourteen known to be part of the nucleosome remodeling deacetylase (NuRD) complex, and four of the Argonaute subfamily) unveiled that they are active. Finally, our results evidenced for the first time a krüppel-associated box (KRAB)-like domain specific to actinopterygians that, together with TRIM33, might allow the functioning of NuRD complex also in fish species. The possible interaction between these two proteins was supported by structural prediction analyses.

Morphological analysis of erythrocytes of an Antarctic teleost under heat stress: Bias of the stabling effect.

The stenothermal Antarctic fish that live in the coastal waters of the Terra Nova Bay (Ross Sea) are rarely exposed to temperatures above zero during the year. We tested whether a slight temperature rise of 1.5 °C affects a sensitive biomarker such as erythrocytes morphology in sections of blood pellets of a small demersal notothen. The erythrocytes' shape descriptors showed significant or highly significant differences temporally from the capture of fish to the conclusion of the experiment. Surprisingly, the erythrocyte's morphology did not show significant differences between the two experimental conditions, returning similar results in control fish stabled at -0.9 °C and in the fish treated at +0.6 °C, although the values of the shape descriptors were often lower in the latter. This study demonstrates the critical issues of comparative physiology in the study of extremely sensitive organisms, such as the fish of the High Antarctic Zone. Moreover, the stabling effect inside the aquarium facilities appears to significantly obscure the effects of the experimental heat treatment.

Investigation of the activity of transposable elements and genes involved in their silencing in the newt Cynops orientalis, a species with a giant genome.

Caudata is an order of amphibians with great variation in genome size, which can reach enormous dimensions in salamanders. In this work, we analysed the activity of transposable elements (TEs) in the transcriptomes obtained from female and male gonads of the Chinese fire-bellied newt, Cynops orientalis, a species with a genome about 12-fold larger than the human genome. We also compared these data with genomes of two basal sarcopterygians, coelacanth and lungfish. In the newt our findings highlighted a major impact of non-LTR retroelements and a greater total TE activity compared to the lungfish Protopterus annectens, an organism also characterized by a giant genome. This difference in TE activity might be due to the presence of young copies in newt in agreement also with the increase in the genome size, an event that occurred independently and later than lungfish. Moreover, the activity of 33 target genes encoding proteins involved in the TE host silencing mechanisms, such as Ago/Piwi and NuRD complex, was evaluated and compared between the three species analysed. These data revealed high transcriptional levels of the target genes in both newt and lungfish and confirmed the activity of NuRD complex genes in adults. Finally, phylogenetic analyses performed on PRDM9 and TRIM28 allowed increasing knowledge about the evolution of these two key genes of the NuRD complex silencing mechanism in vertebrates. Our results confirmed that the gigantism of the newt genomes may be attributed to the activity and accumulation of TEs.

Orally administered nano-polystyrene caused vitellogenin alteration and oxidative stress in the red swamp crayfish (Procambarus clarkii).

Nanoplastics (≤100 nm) represent the smallest fraction of plastic litter and may result in the aquatic environment as degradation products of larger plastic material. To date, few studies focused on the interactions of micro- and nanoplastics with freshwater Decapoda. The red swamp crayfish (Procambarus clarkii, Girard, 1852) is an invasive species able to tolerate highly perturbed environments. As a benthic opportunistic feeder, this species may be susceptible to plastic ingestion. In this study, adult P. clarkii, at intermolt stage, were exposed to 100 µg of 100 nm carboxylated polystyrene nanoparticles (PS NPs) through diet in a 72 h acute toxicity test. An integrated approach was conceived to assess the biological effects of PS NPs, by analyzing both transcriptomic and physiological responses. Total hemocyte counts, basal and total phenoloxidase activities, glycemia and total protein concentration were investigated in crayfish hemolymph at 0 h, 24 h, 48 h and 72 h from PS NPs administration to evaluate general stress response over time. Differentially expressed genes (DEGs) in the hemocytes and hepatopancreas were analyzed to ascertain the response of crayfish to PS NP challenge after 72 h. At a physiological level, crayfish were able to compensate for the induced stress, not exceeding generic stress thresholds. The RNA-Sequencing analysis revealed the altered expression of few genes involved in immune response, oxidative stress, gene transcription and translation, protein degradation, lipid metabolism, oxygen demand, and reproduction after PS NPs exposure. This study suggests that a low concentration of PS NPs may induce mild stress in crayfish, and sheds light on molecular pathways possibly involved in nanoplastic toxicity.

The evolutionary puzzle solution for the origins of the partial loss of the Cτ2 exon in notothenioid fishes.

Cryonotothenioidea is the main group of fishes that thrive in the extremely cold Antarctic environment, thanks to the acquisition of peculiar morphological, physiological and molecular adaptations. We have previously disclosed that IgM, the main immunoglobulin isotype in teleosts, display typical cold-adapted features. Recently, we have analyzed the gene encoding the heavy chain constant region (CH) of the IgT isotype from the Antarctic teleost Trematomus bernacchii (family Nototheniidae), characterized by the near-complete deletion of the CH2 domain. Here, we aimed to track the loss of the CH2 domain along notothenioid phylogeny and to identify its ancestral origins. To this end, we obtained the IgT gene sequences from several species belonging to the Antarctic families Nototheniidae, Bathydraconidae and Artedidraconidae. All species display a CH2 remnant of variable size, encoded by a short Cτ2 exon, which retains functional splicing sites and therefore is included in the mature transcript. We also considered representative species from the three non-Antarctic families: Eleginopsioidea (Eleginops maclovinus), Pseudaphritioidea (Pseudaphritis urvillii) and Bovichtidae (Bovichtus diacanthus and Cottoperca gobio). Even though only E. maclovinus, the sister taxa of Cryonotothenioidea, shared the partial loss of Cτ2, the other non-Antarctic notothenioid species displayed early molecular signatures of this event. These results shed light on the evolutionary path that underlies the origins of this remarkable gene structural modification.

Cold Adaptation in Antarctic Notothenioids: Comparative Transcriptomics Reveals Novel Insights in the Peculiar Role of Gills and Highlights Signatures of Cobalamin Deficiency.

Far from being devoid of life, Antarctic waters are home to Cryonotothenioidea, which represent one of the fascinating cases of evolutionary adaptation to extreme environmental conditions in vertebrates. Thanks to a series of unique morphological and physiological peculiarities, which include the paradigmatic case of loss of hemoglobin in the family Channichthyidae, these fish survive and thrive at sub-zero temperatures. While some of the distinctive features of such adaptations have been known for decades, our knowledge of their genetic and molecular bases is still limited. We generated a reference de novo assembly of the icefish Chionodraco hamatus transcriptome and used this resource for a large-scale comparative analysis among five red-blooded Cryonotothenioidea, the sub-Antarctic notothenioid Eleginops maclovinus and seven temperate teleost species. Our investigation targeted the gills, a tissue of primary importance for gaseous exchange, osmoregulation, ammonia excretion, and its role in fish immunity. One hundred and twenty genes were identified as significantly up-regulated in Antarctic species and surprisingly shared by red- and white-blooded notothenioids, unveiling several previously unreported molecular players that might have contributed to the evolutionary success of Cryonotothenioidea in Antarctica. In particular, we detected cobalamin deficiency signatures and discussed the possible biological implications of this condition concerning hematological alterations and the heavy parasitic loads typically observed in all Cryonotothenioidea.

Transposable Elements and Teleost Migratory Behaviour.

Transposable elements (TEs) represent a considerable fraction of eukaryotic genomes, thereby contributing to genome size, chromosomal rearrangements, and to the generation of new coding genes or regulatory elements. An increasing number of works have reported a link between the genomic abundance of TEs and the adaptation to specific environmental conditions. Diadromy represents a fascinating feature of fish, protagonists of migratory routes between marine and freshwater for reproduction. In this work, we investigated the genomes of 24 fish species, including 15 teleosts with a migratory behaviour. The expected higher relative abundance of DNA transposons in ray-finned fish compared with the other fish groups was not confirmed by the analysis of the dataset considered. The relative contribution of different TE types in migratory ray-finned species did not show clear differences between oceanodromous and potamodromous fish. On the contrary, a remarkable relationship between migratory behaviour and the quantitative difference reported for short interspersed nuclear (retro)elements (SINEs) emerged from the comparison between anadromous and catadromous species, independently from their phylogenetic position. This aspect is likely due to the substantial environmental changes faced by diadromous species during their migratory routes.

Massive gene presence-absence variation shapes an open pan-genome in the Mediterranean mussel.

BACKGROUND: The Mediterranean mussel Mytilus galloprovincialis is an ecologically and economically relevant edible marine bivalve, highly invasive and resilient to biotic and abiotic stressors causing recurrent massive mortalities in other bivalves. Although these traits have been recently linked with the maintenance of a high genetic variation within natural populations, the factors underlying the evolutionary success of this species remain unclear. RESULTS: Here, after the assembly of a 1.28-Gb reference genome and the resequencing of 14 individuals from two independent populations, we reveal a complex pan-genomic architecture in M. galloprovincialis, with a core set of 45,000 genes plus a strikingly high number of dispensable genes (20,000) subject to presence-absence variation, which may be entirely missing in several individuals. We show that dispensable genes are associated with hemizygous genomic regions affected by structural variants, which overall account for nearly 580 Mb of DNA sequence not included in the reference genome assembly. As such, this is the first study to report the widespread occurrence of gene presence-absence variation at a whole-genome scale in the animal kingdom. CONCLUSIONS: Dispensable genes usually belong to young and recently expanded gene families enriched in survival functions, which might be the key to explain the resilience and invasiveness of this species. This unique pan-genome architecture is characterized by dispensable genes in accessory genomic regions that exceed by orders of magnitude those observed in other metazoans, including humans, and closely mirror the open pan-genomes found in prokaryotes and in a few non-metazoan eukaryotes.

Differential Effects of Dietary Supplementation of Krill Meal, Soybean Meal, Butyrate, and Bactocell® on the Gene Expression of Atlantic Salmon Head Kidney.

The head kidney is a key organ that plays a fundamental role in the regulation of the fish immune response and in the maintenance of endocrine homeostasis. Previous studies indicate that the supplementation of exogenous dietary components, such as krill meal (KM), soybean meal (SM), Bactocell® (BA), and butyrate (BU), can have a significant effect on the immune function of the head kidney. The aim of this study was to investigate the differential effect of these four dietary ingredients on the transcriptional profiles of the head kidney of the Atlantic salmon. This study revealed that just a small number of genes were responsive to the feeding regime after a long-term (12 weeks) treatment, and evidenced that the most significant alterations, both in terms of the number of affected genes and magnitude of changes in gene expression, were detectable in the BU- and KM-fed groups compared with controls, while the SM diet had a nearly negligible effect, and BA had no significant effects at all. Most of the differentially expressed genes were involved in the immune response and, in line with data previously obtained from pyloric caeca, major components of the complement system were significantly affected. These alterations were accompanied by an increase in the density of melanomacrophage centers in the KM- and SM-fed group and their reduction in the BU-fed group. While three types of dietary supplements (BU, KM, and SM) were able to produce a significant modulation of some molecular players of the immune system, the butyrate-rich diet was revealed as the one with the most relevant immune-stimulating properties in the head kidney. These preliminary results suggest that further investigations should be aimed towards the elucidation of the potential beneficial effects of butyrate and krill meal supplementation on farmed salmon health and growth performance.

PLANiTS: a curated sequence reference dataset for plant ITS DNA metabarcoding.

DNA metabarcoding combines DNA barcoding with high-throughput sequencing to identify different taxa within environmental communities. The ITS has already been proposed and widely used as universal barcode marker for plants, but a comprehensive, updated and accurate reference dataset of plant ITS sequences has not been available so far. Here, we constructed reference datasets of Viridiplantae ITS1, ITS2 and entire ITS sequences including both Chlorophyta and Streptophyta. The sequences were retrieved from NCBI, and the ITS region was extracted. The sequences underwent identity check to remove misidentified records and were clustered at 99% identity to reduce redundancy and computational effort. For this step, we developed a script called 'better clustering for QIIME' (bc4q) to ensure that the representative sequences are chosen according to the composition of the cluster at a different taxonomic level. The three datasets obtained with the bc4q script are PLANiTS1 (100 224 sequences), PLANiTS2 (96 771 sequences) and PLANiTS (97 550 sequences), and all are pre-formatted for QIIME, being this the most used bioinformatic pipeline for metabarcoding analysis. Being curated and updated reference databases, PLANiTS1, PLANiTS2 and PLANiTS are proposed as a reliable, pivotal first step for a general standardization of plant DNA metabarcoding studies. The bc4q script is presented as a new tool useful in each research dealing with sequences clustering. Database URL: https://github.com/apallavicini/bc4q; https://github.com/apallavicini/PLANiTS.

Molecular Diversity of Mytilin-Like Defense Peptides in Mytilidae (Mollusca, Bivalvia).

The CS-αß architecture is a structural scaffold shared by a high number of small, cationic, cysteine-rich defense peptides, found in nearly all the major branches of the tree of life. Although several CS-αß peptides involved in innate immune response have been described so far in bivalve mollusks, a clear-cut definition of their molecular diversity is still lacking, leaving the evolutionary relationship among defensins, mytilins, myticins and other structurally similar antimicrobial peptides still unclear. In this study, we performed a comprehensive bioinformatic screening of the genomes and transcriptomes available for marine mussels (Mytilida), redefining the distribution of mytilin-like CS-αß peptides, which in spite of limited primary sequence similarity maintain in all cases a well-conserved backbone, stabilized by four disulfide bonds. Variations in the size of the alpha-helix and the two antiparallel beta strand region, as well as the positioning of the cysteine residues involved in the formation of the C1-C5 disulfide bond might allow a certain degree of structural flexibility, whose functional implications remain to be investigated. The identification of mytilins in Trichomya and Perna spp. revealed that many additional CS-αß AMPs remain to be formally described and functionally characterized in Mytilidae, and suggest that a more robust scheme should be used for the future classification of such peptides with respect with their evolutionary origin.

Lysozyme-Induced Transcriptional Regulation of TNF-α Pathway Genes in Cells of the Monocyte Lineage.

Lysozyme is one of the most important anti-bacterial effectors in the innate immune system of animals. Besides its direct antibacterial enzymatic activity, lysozyme displays other biological properties, pointing toward a significant anti-inflammatory effect, many aspects of which are still elusive. Here we investigate the perturbation of gene expression profiles induced by lysozyme in a monocyte cell line in vitro considering a perspective as broad as the whole transcriptome profiling. The results of the RNA-seq experiment show that lysozyme induces transcriptional modulation of the TNF-α/IL-1ß pathway genes in U937 monocytes. The analysis of transcriptomic profiles with IPA® identified a simple but robust molecular network of genes, in which the regulation trends are fully consistent with the anti-inflammatory activity of lysozyme. This study provides the first evidence in support of the anti-inflammatory action of lysozyme on the basis of transcriptomic regulation data resulting from the broad perspective of a whole-transcriptome profiling. Such important effects can be achieved with the supplementation of relatively low concentrations of lysozyme, for a short time of exposure. These new insights allow the potential of lysozyme in pharmacological applications to be better exploited.